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MAGNETIC RESONANCE FLUID ANALYSIS 
APPARATUS AND METHOD 



This application claims priority of provisional application Ser. No. 60/279,324, filed 
5 March 28, 2001 . The content of the above application is incorporated herein by reference. 



FIELD OF THE INVENTION 

The present invention relates gmerally to downhole fluid analysis and more 
particularly to a nuclear magnetic resonance (NMR) apparatus and method providing in-situ 
1 0 data about formation fluids at true reservoir conditions. 

BACKGROUND OF THE INVENTION 

The analysis of downhole fluid samples is desirable in many oil mdustry 

applications. In the prior art this it is typically done by bringing samples to the surface 
15 using sealed containers, and sending the samples for laboratoiy measurements. A number 

of technical and practical limitations are associated with this approach. 

The main concern usually is that the sample(s) taken to the surface may not be 

representative of the downhole geologic formation due to the fact that only limited sample 

material from a limited number of downhole locations can be extracted and taken to the 
20 surface. Thus, taking samples to the surface is impractical if it is desired to measure the 

fluid on a dense grid of sample points. Therefore, by necessity the measurements will only 

provide an incomplete picture of the downhole conditions. 

In addition, these samples firequently contain highly flammable hydrocarbon 

mixtures under pressure. Depressurizing the containers frequently leads to the loss of the 
25 gas content. Handlingof such test samples can be hazardous and costly. Significant 

practical problems are further caused by fluid phase changes during retrieval and transport, 

difficulties in re-creating reservoir conditions and, in general, the significant time delay 

associated with the laboratory analysis. 

It is apparent that the above fluid sampling and analysis process ideally would be 
30 preceded or even completely replaced by downhole fluid analysis of as many samples as 

desired, with the final results instantaneously available at the well site. Nuclear magnetic 

resonance (NMR) technology is well suited for this puipose, as it enables the user to 

determine many properties of an in^situ formation fluid without extracting numerous 

samples. These properties include hydrogen density, self-difiusivity, and relaxation times, 
35 T, and Tj. NMR devices, methods and pulse sequences for use in logging tools are 

described, for example, in U.S. Pat. 4,350,955; 5,557,201; 4,710,713; 4,717,876; 4,717,877; 
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4,717,878; 5,212,447; 5,280,243; 5,309,098; 5,412,320; 5,517,115, 5,557.200; 5,696,448; 
5,936,405; 6,005,389;6,023,164 and 6,051,973. The above patents are hereby incorporated 
by reference. 

5 Direct downhole measurement of certain fluid properties is known in the art. 

Several commercially available tools can be used to this end. Examples include the RDT™ 
tool manufactured by Halliburton, the Reservoir Characterization Instrument (RCf™) from 
Western Atlas, and the Modular Formation Dynamics Tester (MDT™) made by 
Schlumberger. These tester tools have modular design that allows them to be reconfigured 

10 at the well site. Typically, these tools provide pressure-volume measurements, which can 
be used to differentiate liquids from gases, and are also capable of providmg temperature, 
resistivity and other mechanical or electrical measurements. However, these tools do not 
generally provide NMR measurements, as discussed above. 

The use of NMR measurements to determine downhole formation fluid properties is 

15 also known in the field. The first approach to this end is disclosed in U.S. Pat. No. 
6,1 1 1,408 to one of the inventors of the present invention, which discloses a method and 
apparatus for making direct downhole >JMR measurements of formation fluids. The 
disclosure of this patent is incorporated by reference for all purposes. The device m the 
'408 patent, however, generally requires that a portion of the fluid be diverted from the 

20 main flow line and be held stationary for the duration of the measurement, which may take 
about a minute. A possible concern about such use of this device is the occurrence of fluid 
phase separation due to the diversion from the main flow line and separation due to gravity 
once the sample has been stagnated. In addition, depending on relative concentrations, the 
sample chamber may contain only a subset of the phases contained in the flow line. 

25 Accordingly, it is perceived that there is a need for a tester capable of performing 

direct, continuous-flow downhole NMR measurements that can be used to enhance the 
quality and reliability of formation evaluation obtained using prior art techniques. 
Additionally, there is a need to provide a modular NMR downhole analysis apparatus that 
can be used as an add-on to existing testing eqvdpment so as to minimiz e the cost of the 

30 extra measurements. 

In order to more fully appreciate why NMR is important in measurements of this 
type, the reader is directed to Stokes* equation: 

DockT/Ti, (k=1.38xlO^^J/K) (1) 

35 
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Eq. (1) essentially indicates that the self-difiusion coefficient D is mversely 
proportional to viscosity i] and vice versa. Viscosity and diflEiisivity are both related to the 
translational motion of molecules and therefore must be interrelated. At higher 
temperatures T, a molecule contains more energy and can move fester against a given 

5 "friction'* r\, therefore D oc T. Diffusivity is a property that can be precisely determined by 
NMR techniques without disturbing or altering the fluid. The relationship D oc T/r\ has 
. been verified over a wide range of viscosities at different temperatures and pressures by 
NMR spin-echo experiments. Referring to the prior art references listed below, see, for 
example, Abragam, 1 96 1 . 

1 0 One has to be more careful with relationships involving the NMR relaxation times 

Tj and Tj. The applicability of expressions of the form, 

T„TjOckT/Ti, • (2) 

^5 is more limited than that of Eq. (1). The main reason is that gas/liquid mixtures, such as 
live oils, relax by more than one relaxation mechanism: dipole-dipole for the liquid phase 
and mainly spin rotation for the gas phase. In combination, however, NMR relaxometry 
(measuring Tj and Tj) and NMR diffusometry (deteraiining D) are powerful tools to 
characterize live oils. 

20 With further reference to the prior art list below, the study of NMR relaxation times 

with respect to oil properties began shorfly after NMR was first demonstrated (Bloembergen 
et al, 1948; Brown, 1961). The practical aspects of how to relate Ti, Tj and D to 
petrophysical fluid properties such as viscosity and gas/oil ratio have attracted interest much 
more recentty (Klemberg et aU 1996; Lo et a/., 1998; Zhang et a/., 1998; Appel et aL, 

25 2000; Lo et aL, 2000). These investigations are significant because NMR relaxometry can 
be performed at much lower field strengths and with much lower homogeneity requirements 
than NMR spectroscopy (1,000 ppmv.<l ppm). The latter seems to be out of reach for 
downhole applications, but the development of a robust and accurate NMR 
relaxometiy/diffusometiy system for downhole use proved to be feasible. 

30 Another favorable factor was the development of the Reservoir Description Tool of 

Halliburton (RDT; Proett et al, 1999), which is a modular wireline samplmg and testing 
system that can be configured in a variety of tool combinations and can readily accept add- 
on analysis devices. Advantageously, these tools generally can operate independently of 
each other. The inclusion of an NMR fluid analyzer in the RDT tool string has the 

3^ following benefits: 
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(1) The level of fluid contamination due to, for example, filtrates from water or oil 
based muds can be assessed continuously by observing the fluid's relaxation time 
distribution. 

(2) Estimates for fluid viscosity and gas/oil ratio (GOR) can be estunated in real 

5 tune. Viscosity can be derived from either T| or D, and GOR from a combination of both 
parameters. 

(3) The NMR analysis takes place at true reservoir conditions, removing ambiguities 
associated with sampling and transport procedures. 

(4) The shared measvirement principles behind wireline NMR tools and the 
1 0 downhole NMR fluid analyzer encourage the development of integrated formation 

evaluation methods. 

(5) Data from the NMR analysis are processed downhole and the results are 
viewable in real time. This information can be useful in tiie selection of further sampling 
points and for the fine-tuning an open-hole logging program. A simple example of this 

1 5 synergy is the verification of the hydrogen density m the invaded zone by in-situ fluid 

testing, a number which feeds directiy into the porosity reading of both wirelme and 

logging-while-drilling (LWD) tools. 

Details of the apparatus and method in accordance with the present invention are 

provided below. The interested reader is directed for additional background information to 
20 the disclosiore of the following references, which are incorporated herein by reference for 

background. For simplicity, in tiie foUowmg disclosure only the first author and date of 

publication are provided. 

1 . Abragam, A,, 1961, Principles of Nuclear Magnetism, Oxford University Press, 
Oxford. 

25 2. Appel, M., Freeman, J.J., Perkins, RB., and van Diik, N.P., 2000, Reservoir Fluid 
Study by Nuclear Magnetic Resonance, Paper HH: SPWLA, presented at the 41*' Annual 
Logging Symposium, Dallas, TX, June 4-7. 

3. Bioembergen, N., Purcell, E. M., and Pound, R V., 1948, Phys. Rev. 73, p. 679. 

4. Brown, R J. S,, 1961 , Proton Relaxation in Crude Oils, Nature 189, no. 4762, p. 
30 387. 

5. Coates, G. R., Xiao, L., and Prammer, M. G., 1999, NMR Logging, Prmciples and 
Applications, Halliburton Energy Services, Houston. 

6. Freedman, R., Sezginer, A., Flaum, M., Matteson, A., Lo, S., and Hirasaki, G. J., 
2000, A New NMR Method of Fluid Characterization in Reservoir Rocks: Experimental 

35 
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Confinnation and Simvdation Results, Paper SPE-63214: Society of Petroleum Engineers, 
presented at the 75* Annual Technical Conference and ExWbition, Dallas, TX, Oct 1-4. 
7. Kimmich, R., 1997, NMR Tomogrqphy, Diffiisometry. Relaxometry, Springer- 
Verlag, Berlin. 

5 8, Kleinberg, R. L., and Vinegar, H. J., 1996, NMR Properties of Reservoir Fluids, The 
Log Analyst, Nov.-Dec,, p. 20. 

9. Lo, S.-W., Hirasald, Kobayashl, R., and House, W. V., 1998, Relaxation Time and 
Difiusion Measurements of Methane and n-Decane Mixtures, The Log Analyst, Nov.-Dec, 
p. 43. 

10 10. Lo, S.-W., Hirasaki,G., House, W. v., and Kobayashi,R., 2000, Correlations of 
NMR Relaxation Times with Viscosity, Diffusivity, and Gas/Oil Ratio of 
Methane/Hydrocarbon Mixtures, Paper SPE-63217: Society of Petroleum Engineers, 
presented at the 75* Annual Technical Conference and Exhibition, Dallas, TX, Oct. 1-4. 

1 1 . Prammer, M. G., 1 994, NMR Pore Size Distributions and Permeability at the Well 
15 Site, Paper SPE-28368: Society of Petroleum Engineers, presented at the 69* Annual 

Technical Conference and Exhibition, New Orleans, LA, Sep. 25-28. 

12. Prammer, M. G., Mardon, D., Coates, G. R., and Miller, M. N, 1 995, Lithology- 
Independent Gas Detection by Gradient-NMR Logging, Paper SPE-30562: Society of 
Petroleum Engineers, presented at 70* Annual Technical Conference and Exhibition, 

20 Dallas, TX, Oct. 22-25. 

13. Prammer, M. G., Bouton, J., Chandler, R. N., Drack, E. D., and Miller, M. N., 1998, 
A New Multiband Generation of NMR Logging Tools, Paper SPE-4901 1 : Society of 
Petroleum Engineers, presented at the 73"* Annual Technical Conference and Exhibition, 
New Orleans, LA, Sept. 27-30. 

25 14. Pranuner, M. G., Drack, E., Goodman, G., Masak, P., Menger, S., Morys, M., 
Zannoni, S., Suddarth, B., and Dudley, J., 2000, The Magnetic Resonance While-Drilling 
Tool: Theory and Operation, Paper SPE-62981 : Society of Petroleum Engineers, presaited 
at the 75* Annual Technical Conference and Exhibition, Dallas, TX, Oct. 1-4. 

1 5. Prammer, 2000, Method and apparatus for differentiating oil based mud filtrate from 
30 connate oil, U.S. PaL 6,107,796, Aug. 22. 

1 6. Proett, M. A., Gilbert, G. N., Chin, W. C, and Monroe, M. L., 1999, New Wireline 
Formation Testing Tool with Advanced Sampling Technology, Paper SPE-5671 1, Society 
of Petroleum Engineers, presented at the 74* Annual Technical Conference and Exhibition, 
Houston, TX, Oct. 3-6. 
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17, Van Dusen, A., Williams, S., Fadnes, F. H., and Irvine-Fortescue, J,, 2000, 
Determination of Hydrocarbon Properties by Optical Analysis During Wireline Flxaid 
Sampling, Paper SPE-63252: Society of Petroleum Engineers, presented at the 75* Annual 
Technical Conference and Exhibition, Dallas, TX, Oct. 1-4. 
5 18. Zhang, Q., Lo, S., Huang, C. C, Hirasaki, G. J., Kobayashi, R., and House, W. V., 
1998, Some Exceptions to Default NMR Rock and Fluid Properties, Paper FF: SPWLA, 
presented at the 39* Annual Logging Symposium, Keystone, CO, May 26-29. 

10 SUMMARY OF THE INVENTION 

The NMR analyzer apparatus and method of the present invention provides in-situ 
data of foraiation fluids at true reservoir conditions and overcome the above-identified and 
other problems associated with the prior art. 

In one aspect, the present invention is a method for analyzing formation fluids in a 

15 borehole environment, comprising the steps of: (a) kitroducing formation fluids in a flow- 
through vessel located in the borehole; (b) generatmg a substantially imiform static 
magnetic field in the vessel with a defined magnetic field direction; (c) generating pulsed 
oscillating magnetic fields for exciting nuclei of formation fluids in a first portion of the 
vessel, said oscillating fields having a magnetic direction substantially perpendicular to the 

20 direction of the static magnetic field; (d) receiving nuclear magnetic resonance (NMR) 
relaxation signals fi-om excited nuclei of the fluid in a second portion of the vessel, smaller 
than the first portion; and (e) analyzing the received signals to determine fluid properties at 
any flow rate of the formation fluids in the vessel below a predetermined non-zero 
threshold. In a specific embodiment, the pulsed magnetic fields used are according to a 

25 saturation recovery pulse sequence. In various other specific embodiments, the threshold 
value for the flow rate is determined by the length of the first portion, the second portion of 
the vessel is positioned downstream &om the first portion of the vessel, and the first and 
second portions of the vessel overlap at least in part. In a q)ecific application, the steps of 
the method are performed substantially continuously, and the step of analyzing comprises 

30 assessing mud filtrate contamination for the* formation fluids. Additionally , the step of 
analyzing may comprise monitoring the Ti profile of fluids passed through the vessel. 

In another aspect, the present invention is a method for analyzing fluids, comprising 
the steps of: (a) providing a flow-through passage for fluids m a measurement chamber, the 
fluids having flow rate within a pre-determmed range; (b) performing aNMR experiment to 

35 excite substantially all nuclei of the fluids in the chamber at a given time interval; and (c) 
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processing NMR signals obtained from a portion of the chamber in said experiment without 
regard for the flow rate of the flxiids in the chamber. Various specific embodiments are 
disclosed in detail next. 

In another aspect, the present invention is an apparatus for analyzing downhole 

5 formation fluids, comprising: (a) a conduit for introducing formation fluids into the 
apparatus and for providing flow-through passage, the conduit having an inlet end and an 
outlet end; (b) at least one magnet assembly enclosing the conduit for generating in the 
condmt a substantially uniform static magnetic field with a defined magnetic field direction; 
(c) at least one transmitting antetma operative to generate pulsed magnetic fields in the 

1 0 conduit in a direction substantially perpendicular to the static field direction for exciting 
nuclei of fluids contained in the condmt; and (d) at least one receiving antenna operative to 
receive NMR signals from fluids in a portion of the conduit, the receivmg antenna being 
shorter than the at least one transmitting coil, so that received NMR signals correspond only 
to a portion of the excited nuclei. The apparatus of this aspect may fiorther comprise a 

1 5 shield mounted between the transmitting anteima and a portion of the magnet assembly. In 
various specific embodiments, the apparatus may have a portion of the conduit between the 
inlet and outlet ends be wider in dimension dian the conduit at either end, and finther have a 
diffuser positioned near the inlet end of the conduit, for providing consistent fluid flow 
velocity over the wider portion of the conduit. In a preferred embodiment, the conduit is 

20 adapted for attachment to the flow line of a modular wireline logging tool. Additionally, the 
at least one magnet assembly in a preferred embodiment has a polarization portion located 
near the inlet end, and a resonance portion located near the outlet end, the strength of the 
magnetic field in the polarization portion being higher than the strength of the magnetic 
field in the resonance portion. Other features of the preferred embodiments are disclosed in 

25 more detail next. 

In a fiirther aspect, flie present invention is a tester module for use with modular 
downhole formation testers for downhole NMR testing of formation fluids comprismg: a 
vessel for providing flow-through passage of formation fluids and for conducting downhole 
NMR measurements, said vessel being adapted to withstand borehole environment 

30 conditions; at least one tubular magnet defining a longitudinal axis, the magnet having 
magnet sections with magnetization direction(s) perpendicular to the longitudmal axis, and 
enclosing the passage to generate therein a static magnetic field with predetermined 
magnetic field direction; at least one radio frequency (RF) transmitter operative to generate 
pulsed RF magnetic fields m a first portion of the passage m a direction substantially 

35 perpendicular to the static field direction for exciting nuclei of fluids in the passage; a 
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receiver for acqxiiring NMR signals from excited nuclei in a second portion of the passage 
smaller than the &st portion and positioned do>vnstream therefrom, and a processor for 
analyzing properties of fluids in the passage based on signals from the receiver. 
In yet another aspect, the present invention is a n apparatus for downhole 

5 differentiating between fluid types present in a geologic formation comprising: (a) a 
measurement chamber having flow-through passage for formation fluids; (b) an NMR 
testing module capable of performing an* NMR experiment on a portion of the formation 
fluid witiiin the passage, the experiment exciting nuclei of the formation fluids in tiie 
passage at a given time interval; and (c) processor receiving NMR signals obtained in said 

1 0 testing module to diflferentiate fluid types in the foraxation fluids without regard for the flow 
rate of the formation fluids in the passage. 

Various specific features of the preferred embodiments are disclosed in detail below 
and are defined in the appended claims. 

15 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure lis a schematic diagram of a downhole NMR fluid analyze: in accordance 
with a preferred embodiment of the present invention; Fig. 1 A is a simplified diagram; 
5 Figures 2A-D illustrate horizontal cross-sections of the apparatus shown in Fig. 1. 

Figure 3 illustrates a schematic block diagram of the electronics used in accordance 
with one embodiment of the present invention; 

Figure 4 illustrates a saturation-recovery pulse sequence diagram used for T, 
measurements in accordance with a preferred embodiment of the present invention; 
10 Figure 5 illustrates examples of T, saturation-recovery data for three different fluids 

as seen in the NMR fluid analyzer in accordance with the present invention; 
Figure 6 illustrates the data shown in Fig. 5, but in the T, domain; 
Figure 7 illustrates a pulse sequence employed in accordance vdth a preferred 
embodiment for diffusivity measurements; 
1 5 Figure 8 is an example of a diffusivity measurement in accordance with the present 

invention. 
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pFT^n.T^n nRSri^TPT TON OF THE PREFERRED EMBODIMENTS 

The Apparatus 

Figure 1 shows a schematic diagram of a fluid analysis apparatus 100 used in 

5 accordance vdth a preferred embodiment of the present invention. A simplified diagram of 
the analysis apparatus of the invention is shown in Fig. 1 A. The apparatus of this invention 
is advantageously designed as a module for use with a logging tool, such as for example the 
RDT™ tool string by Hallibiirton Energy Services Corporation, which tool is configurable 
for a variety of sampling, testing and monitoring purposes. 

10 As shown, fluids enter the device at the top and pass through two sections, referred 

to as polarization and resonance sections, respectively. Measurements are performed as the 
fluid flow 1 5 passes through the device. The fluid entering the system is initially subjected 
to a strong magnetic field to achieve rapid polarization of the hydrogen nuclei. NMR 
measurements take place in the lower section, where the field strength is lowered. In a 

1 5 preferred embodiment, two separate radio frequency coils are used for pulse transmission 
and for reception. This split scheme allows for a transmitter coil 30 that is longer than the 
receiver 35. In accordance with this invention, by pulsing a larger volume of fluid and by 
receiving only from the bottom portion, relaxation times can be determined without regard 
for the actual flow rate. 

20 More specifically, with reference to Fig. 1, the fluid to be analyzed flows into the 

measurement chamber 10, which is wddened to a diameter of approximately 2 cm in a 
specific embodiment. The increase in diameter (compared to a standard 0,2 inch - i.e., 

0. 508 cm - flow tube diameter) serves two purposes: it increases the available volume and 
NMR signal amplitude and it slows down the fluid, which allows to shorten the length of 

25 the fluid analyzer apparatus. In a specific embodiment, the incoming fluids encounter a 
conical diffuser 5 that breaks the typical plume-like flow pattern which develops when a 
narrow flow channel widens into a larger volume. The result is a more consistent flow 
velocity over the entire cross section of the measurement volume. The length of the 
analysis apparatus is dictated by the requirement for suflBcient polarization time at various 

30 possible pump-out rates, and in a preferred embodiment is about 1 meter. 

In alternative embodiments the measxirement volume may be filled with collimators, 

1. e., thin tubes, the goal being to slow down the fluid as much as possible across the entire 
cross section, while maintaining the initial flow conditions. 

In a specific embodiment, fluids can be pressurized to up to 25,000 psi, although 
35 higher or lower pressure values can be used in alternative implementations. Under such 
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high pressure conditions, there must be outside counter-pressure, which in the above 
specific embodiment is selected to be about 20,000 psi, so that the structure is designed to 
absorb a diflferential pressure of about 5,000 psi. Most of this differential pressure is 
absorbed in a steel hull (not shown), serving as a pressure barrel that in a preferred 
5 embodiment has connectors on each end that mate to a tool string, for example of the RDT 
tool. Preferably, inside components of the device are either solid (magnetic elements) or are 
potted in solid compounds to allow the pressure to pass through and not to damage any 
components. 

As shown in Fig. 1 , in a preferred embodiment there are four main sections of the 

10 analyzer apparatus, labeled (A) through (D). Sections (A) and (B), which occupy tiie top 
portion of the instrument are required to rapidly polarize and then stabilize the nuclear 
magnetic moment. The equilibrium magnetization of the device in this embodiment 
corresponds to an applied field of 1,000 Gauss, although up to about 2,500 Gauss can be 
applied in section (A) in order to achieve very fast polarization. The NMR measurement 

15 takes place in sections (C) and (D), where section (C) is longer than D. In the specific 
embodiment, section (C) is about two-three times the length of section (D). As noted, the 
radio frequency coil of the device is preferably split into a transmit-only portion (C) and a 
transmit/receive portion (D). In accordance with the mvention, it is important that the fluid 
flows from section (C) into (D) and not vice versa, because otherwise measurements would 

20 be inaccurate. Thus, excitation pulses are applied to tiie entire resonating volvune (C)+(D), 
while reception is only made from section (D). As described in fiirther detail below, the 
nature of a saturation-recovery Ti experiment, used in accordance with a preferred 
embodiment of the method, is that a single pulse can be used to prepare a large volume, 
while the actual readout may happen over a smaller volume. In the present invention, this 

25 feature is used to render the measurement of T^ relaxation times largely independent of flow 
velocities, up to a certain practical limit. The RF antennas in a preferred embodiment are 
wound bifilar, where tiie transmit coil is longer, extending throughout regions (C) and (D). 
In a specific embodiment, transmission and reception both operate at 4.258 MHZ at room 
temperature, consistent with 1,000 Gauss field strengtii. At higher temperatures, the 

30 operating frequency is reduced to track tiie reversible reduction in magnetic field strength. 

As fiuther shown in Fig. 1, the analysis apparatus may also have a permeable flux 
return and magnetic shield 20, magnets in configurations discussed in detail below, Faraday 
shields 40 and 45 and oil batii 50. Preferably, in close proximity can be placed one or more 
temperature sensors (not shown) that track the magnets' temperature for frequency 

35 adjustment and the sample's temperature for viscosity calculations. As described next. 
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using the apparatus illustrated in Fig. 1 it is possible to determine hydrogen density, 
relaxation times and self-diffusivity without altering the fluid. 

The magnetic field in the measurement volume of tiie device shown in Fig. 1 in 
general is not entirely uniform. This volume can be split conceptually into an interior 
5 region, where the field gradient is negligible, and a fiinge region, where the field changes 
with an approximately uniform gradient. The fiinge region may comprise about one third of 
the total sensitive volume. During T, measurements and at short pulse-to-pulse spacings 
(such as 0.25 ms), the effect of the gradient is not noticeable. To perform a diffusion 
measurement, m accordance with a preferred embodiment the main flow is diverted and a 
10 sample is stagnated within the NMR chamber. Furthermore, the pulse-to-pulse spacing (Te) 
is increased to induce diffusion-dependent signal dephasing. The imiform and the fiinge 
regions are large compared to the largest possible diffusion length, therefore, these regions 
are essentially isolated from each other for the -duration of a single pulse-echo train. 
Fig. 2A is a cross-sectional view of the first section A of the apparatus that the 
15 flowing fluid encounters after it enters the measurement chamber. The purpose of this pre- 
polarization section is to polarize hydrogen nuclei in the fluid(s) as rapidly as possible, so 
that they will exhibit fiill polarization under the operating field. The flow tube 60, is made 
from ceramic, glass or PEEK material, and is surrounded by a Faraday shield 45. This 
shield consists in a preferred embodiment of a copper cylinder tube, slotted vertically along 
20 its axis. In a preferred embodiment, the magnet consists of eight segments 25 that are 
magnetized as indicated by the anows in Fig. 2A. The magnet material is preferably 
Samarium-Cobalt with very low temperature coefficient. A supplier for this material is 
Electron Energy Corp., Lancaster, Pa. This magnet configuration is known as the 
"Halbach" configuration, which produces a highly uniform and strong field in the interior. 
25 The magnet is dimensioned to produce a field of about 2,000-2,500 Gauss within the flow 
tube 60. The exact value is not critical and may vaiy by several hundred Gauss. This 
magnet is encased in a magnetically permeabjie steel tube that acts as magnetic shield and 
confines the field to the inside of the apparatus. 

Fig 2B is a cross-sectional view of the second section (B) of the apparatus, which is 
30 located between the polarization section (A) and the measurement sections (C) and (D). 
The purpose of section (B) is to allow the hydrogen spins to setUe to an equilibrium 
polarization that is close to a non-flowing magnetization correspondmg to an external field 
of 1,000 Gauss. Without this section it would be possible that overpolarized fluids may 
enter the measurement section below and cause a distortion in the NMR amplitude tiiat is 
35 dependent on the flow rate. In a preferred embodiment, the transversal magnetic field 
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generated in the flow tube is approximately 850 Gauss, but again the exact value is not 
critical and can be varied in practical implementations. As shown, all magnet elements of 
this section are magnetized in the same direction, resulting in a weaker field than that in 
section (A). 

5 Fig. 2C is a cross-sectional view of the third section (C) of the apparatus. This is 

the transmit-only section that follows the stabilization section (B). In a preferred 
embodiment, the transversal magnetic field is about 1,000 Gauss, accurate to +/- 10 Gauss, 
corresponding to a hydrogen NMR resonance fi:equency of 4.26 MHz. A solenoidal coil 30 
woimd around the measurement volume generates an RF field of fliat firequency in 

10 longitudinal du-ection (out of the paper plane). The solenoid is electrically shielded against 
the fluid by a slotted copper tube 45. This tube terminates any electric field emanating fi:om 
the coil. The shield allows magnetic flux to pass through as long as there is at least one slot 
running jfrom top to bottom parallel to the longitudinal axis. A second shield in the form of 
a second copper tube 40 is mounted between the solenoid and the magnet elements. This 

1 5 shield is not slotted, because it needs to pass only a static magnetic field. All dynamic 
electromagnetic fields are terminated at this shield as well as any interference or noise firom 
the outside. Any fluid(s) within this section typically are fiilly polarized corresponding to 
the working field of 1,000 Gauss. The application of a firequency-swept pulse that covers 
the frequency range corresponding to fields fi:om 990 to 1010 Gauss saturates the nuclear 

20 magnetization. In accordance with the present invention the subsequent built-up towards 
the equilibriimi condition can be monitored in section (D) of the apparatus. 

Fig, 2D is a cross-sectional view of the last section (D) of the apparatus 100. This 
is the transmit/receive section of the device. In a preferred embodiment, the transversal 
magnetic field is 1,000 Gauss, accurate to +/- 1 Gauss, corresponding to a hydrogen NMR 

25 resonance fi-equency of 4.26 MHz. A solenoidal coil 35 wound around the measurement 
volume generates an RF field of that frequency in longitudinal direction (out of the paper 
plane) and also receives induced NMR signals from the measurement volume. The two 
Faraday shields are identical to those in section (C). The NMR time constants T, of the 
fluid(s) under investigation are determined by varying the delay time between a broadband 

30 saturation pulse and a read-out sequence. If the flow velocity does not exceed 10 cm/s, the 
measurement is generally independent of the flow speed and of the flow profile. 

In manufacturing this magnet system, a large number of segments is produced, 
magnetized and their actual magnetic remanence Br is measured. The best-matching 
segments are preferably used to build up section (D), such that the field variations are as 

35 small as possible, preferably not to exceed 0.1%. Segments with inferior match are used for 
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section (C), where the accuracy requirement is reduced to about 1%. The magnetic 
elements for sections (A) and (B) have different dimensions and are uncritical in their 
magnetic remanence. 

The electronics used in the analysis apparatus of the present invention is similar to 

5 that of a NMR spectrometer and is illustrated in a block diagram form in Fig. 3. The 
comparatively low frequency of 4.2 MHz used in a preferred embodiment of the invention 
allows many traditionally analog functions to be realized readily as digital signal processing 
(DSP) algorithms. A frequency source 34, controlled by a pulse programmer, sends its 
signal to a power amplifier 37, which in turn drives the transmitter antenna 30. All timing 

10 functions, like pulse widths and acquisition windows, are ftiUy prograinmable. On the 
receive side, the signal from the receiver antenna 35 is amplified, synchronously 
demodulated and integrated. In a preferred embodiment the system also performs its own 
calibration. All pertinent calibration factors are stored in the tool itself and after calibration 
echo amplitudes are reported on a scale of 0-2. 

15 More specifically, as shown in Fig. 3, the two coils of the device 30, 35 are 

connected to resonating capacitors 31 . These capacitors are of the NPO (no temperature 
coefficient) and PTC type (positive temperature coefficient), shunted in parallel, as shown. 
The resultant temperature characteristic is such that with increasing temperature, when the 
static magnetic field weakens (typically 1% per 100°C), the capacitance increases at twice 

20 the rate (typically 2% per 100°C). The resultant LC circuit resonant frequency drops at half 
the capacitor rate (1% per lOOT) and therefore follows the NMR resonance, making re- 
tuning of the circuit uimecessary. 

In a transmit mode, the controller 33 gates the signal generator 34 of the apparatus 
and the power amplifier 37 to produce a radio frequency pulse in both coils. The high 

25 voltage applied causes all crossed diodes 39 to conduct, thereby connecting the two coils. 
In receive mode, the crossed diodes stop conducting and signal is only received from the 
lower coil 35. The signal is amplified, digitized and fed into the digital signal processor 33 
for demodulation and fiirdier processing. 

30 



35 
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Measurements 

A. Understanding NMR Relaxation Times 

NMR measurements are based on the observation that when an assembly of 

5 magnetic moments, such as those of hydrogen nuclei, are exposed to a static magnetic field 
they tend to align along the direction of the magnetic field, resulting in bulk magnetization. 
The rate at which equilibriimi is established in such bulk magnetization upon provision of a 
static magnetic field is characterized by the parameter Tj, known as the spin-lattice 
relaxation time. The spin-lattice relaxation time Ti describes the coupling of nuclear spins 

10 to energy-absorbing molecular motions like rotation, vibration and translation. In most 
fluids at or above ambient temperatures, the coupling to these modes is very inefficient, 
resulting in T,'s in the millisecond and second range. 

Another related and fiequentiy used NMR parameter is the spin-spin relaxation time 
constant Tj (also known as transverse relaxation time), which is an expression of the 

15 relaxation due to non-homogeneities in the local magnetic field over the sensing volume of 
the logging tool. The mechanisms for spin-spin relaxation time T^ include, in addition to 
those contributing to T„ the exchange of energy between spins. These effects are small in 
buUc fluids and therefore Tj basically equals T^. Spin-spin coupling is relevant in heavy oil 
components, such as asphaltenes, resins, etc. Both relaxation times provide indirect 

20 information about the formation porosity, the composition and quantity of the formation 
fluid, and others. 

Another measurement parameter used in NMR is the formation diffiisivity. 
Generally, difiusion refers to the motion of atoms in a gaseous or liquid state due to their 
thermal energy. Self-diffusion of a fluid is directly related to the viscosity of the fluid, a 

25 parameter ofconsiderable importance in borehole surveys, hiaunifomi magnetic field, 
diflfijsion has littie effect on the decay rate of the measured NMR echoes. In a gradient 
magnetic field, however, diffusion causes atoms to move fi^om their original positions to 
new ones, which also causes these atoms to acquire different phase shifts compared to 
atoms that did not move. This contributes to a faster rate of relaxation. 

30 Expressions for Tj in closed form have been derived only for spherical molecules. 

This is obviously a poor approximation to chain-type hydrocarbon molecules, but should 
serve as illustration for the underlying mechanisms. Translational and rotational relaxation 
rates are, respectively: 



35 1 / T, (translational) « Nti / kT, (3) 
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1 / T, (rotational) « (a^ / b*)Ti / kT, (4) 



where N is the spin density of the fluid (proportional to the hydrogen index Ih), a is the 

molecular radius and b is the distance between hydrogen spins on the same molecule. Both 
S mechanisms have the same basic relationship with viscosity and temperature, which 

explains the relative simplicity of Eq. (2). 

Pressure exerts its effect by changing fluid density. Pure liquids are barely 

compressible and the effect pressure has on them is limited. Downhole pressures over 

10,000 psi can induce significant changes in fluid viscosity (Jones, 1 991). As a rule-of- 
10 thxmib, high pressures and high temperatures have opposing effects. Higher temperatures 

expand fluids and increase mobility and relaxation times, while increased pressures reduce 

mobility and relaxation times in liquids. 

Probably the biggest effect downhole high pressures have is the increase in dissolved 

gas volume. The relaxation times in the gas phase exhibit the following behavior: 

15 

Ti = Tj ocT]/kT (methane vapor) (5) 

The reason the relationship in Eq, (2) is not followed is that protons in the gas phase 
relax by spin-rotation and not by dipole-dipole interaction. Increasing the pressure on gas 
20 increases T„ which is contrary to the behavior of fluids. For the gas signal to be detectable 
downhole, the bulk density and the hydrogen density must be relatively high, i.e. only the 
high- viscosity end of Eq. (5) is visible. Under these circumstances, methane relaxation 
times are in the range from hundreds of milliseconds to several seconds. 

Dissolved gas has a profound effect on the Ti of oil. See, Appel et ai, 2000. Lo et 
25 aL (2000) have developed a mixing-rule model for methane-alkane mixtures that links GOR 
to T, and difiusivity. In this model, T, consists of two components: one proportional to 
kT/Tj (Eq. 2) and one proportional to VkT (Eq. 5). The observed Tj is a combination of the 
two, weighted by the proton fractions of the alkane and the methane gas. 

Based on the disciission above, following are examples of measurements that can be 
30 made in accordance with the method of the present invention. 

B, Relaxation time measurements 

In the apparatus of the present invention Tj relaxation times and hydrogen density 
can be measured continuously whether or not the fluid is stagnant or flowing. This mode is 
35 most useful during the pump-out period to get an initial assessment of mud filtrate 
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contamination. Fig. 4 illustrates the pulse sequence employed used in a prefeired 
embodiment It will be appreciated that it is a standard saturation-recovery sequence, where 
an initial saturation pulse is followed by a variable delay. In a preferred embodiment, the 
delay is programmable and is typically stepped through the values 1, 2, 4, 8, 16, ... , 16384 

5 ms in cyclical £ashion. The recovered magnetization at the end of the delay is determined 
by a short read-out sequence, consisting of two pulses and one spin echo. The height of the 
echo, if plotted as function of delay time, traces out a recovery curve that is converted into a 
T, distribution by standard inversion methods. In one embodiment, the inversion algorithm 
is a variant of the method employed to calculate T2 distributions from wireline data, as 

10 disclosed m U.S. Pat. No. 5,517,1 15. See also Prammer, 1994 reference. With the above 
sequence it takes 33 seconds to complete a measurement cycle. The signal-to-noise figure 
of the system is so high that additional averaging may be unnecessary. The sequence 
described above is insensitive to fluid flow and can be used to continuously monitor the Tj 
profile of pumped fluids. 

15 More specifically, the Tj measurement sequence is initiated by a fi'equency-swept 

saturation pulse. With reference to Fig, 1, its frequency is selected such that the entire range 
of resonance frequencies within the sections (C) and (D) of the apparatus are affected. In a 
specific embodiment, this range is typically the NMR center frequency +/- 1% (4.2 MHz +/- 
40 kHz). Pulse amplitude, frequency sweep rate and pulse length are adjusted to effect 

20 saturation in volumes (C) and (D). 

As noted, following the saturation pulse, a variable delay is inserted. In a preferred 
embodiment, consecutive measurements with delay values of 1 ms, 2 ms, 4 ms, ... up to 
16384 ms are used. During these intervals, the nuclear magnetization builds up again to its 
equilibrium value. Also during this time, depending on the flow rate, fluid volume moves 

25 from volume (C) into volume (D), while unprepared fluid enters volume (C). As long as the 
flow rate is not high enough to allow unprepared fluid from (A) or (B) to enter (D), it will 
be appreciated that the measurement is independent of the actual flow rate. After the 
saturation-recovery delay, the instantaneous value of the nuclear magnetization is 
determined. This is done wdth a short pulse sequence, consisting of a nil pulse, followed by 

30 an pulse. The RF phase of these pulses is shifted by 90° against each other to cancel the 
effects of Bo and RF field imperfections. This is equivalent to the start of a CPMG 
sequence. The time between these pulses is typically 0.125 ms and 0.125 ms after the n 
pulse a spin echo forms. This echo is digitized, quantified and its amplitude is taken as a 
measure of the recovered magnetization as function of the saturation-recovery delay. Note 

35 that the Tc^ and 7C pulses can be narrow-band and need not be frequency-swept. The reason 
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is that they are only relevant for section (D) which has a very tightly controlled field and 
resonance frequency distributioa 

Examples for T| distributions for some example fluids obtained in accordance wiht 
the method of this invention are shown in Figs. 5 and 6. The data points were acquired 

5 according to the sequence in Fig. 4 and inverted from the time domain to the T, domain. 
The horizontal axis, **time " in Fig. 5 denotes the time elapsed between the saturation pulse 
and the readout sequence, the vertical axis is signal amplitude in arbitrary units. The results 
are easier to interpret after inversion from time domain to Tl domain, as shown in Fig. 6. In 
this example 53 points are specified for the inversion result. The single, sharp peak at 2-3 s 

10 is characteristic of water, the rounded peak in the "oil window*' 0.5-1 s indicates oil and the 
broad response from the crude oil in the bottom panel is characteristic for complex 
hydrocarbons. Shown in this figure here are examples of Tj saturation-recovery data for 
three different fluids: brine, Diesel oil and a crude oil. 

The data points illustrated have been acquired by circulating different fluids through 

15 the analyzer. Shown from top to bottom are; water (mild brine) with a single relaxation 
peak in the ' Vater window" at 2 seconds; next a simple hydrocarbon (diesel) with a single 
relaxation peak in the "oil window" at 0.5-1 second; and a complex hydrocarbon (crude), 
which shows a characteristic asymmetric distribution that starts in the few tens of 
milliseconds and extends to the "oil window." These samples were \mder atmospheric 

20 conditions at ambient temperature. At elevated temperatures, Eq. (2) predicts an increase in 
Tl proportional to the absolute temperature in addition to increases due to reduction in 
viscosity. 

It is important to note that using the apparatus and method of the present invention 
the determination of long relaxation times does no longer depend on how long an echo train 
25 persists. In the implementation discussed above small perturbations in the applied field 
have relatively limited effect. Additionally, the saturation pulse prepares a much larger 
sample volume than what is actually used for the readout portion. Therefore, as long as the 
flow rate is low enough, and the readout is based on a fluid sample that was present 
anywhere within the regions (C) or (D) during the saturation pulse, the measurement is 



In contrast to Tj, the Tj parameter generally caimot be determined on a flowing 
sample. Distributions of T2 times are detennined in accordance with the present invention 
by standard Carr-Purcell-Meiboom-GiU (CPMG) sequences on samples that have been 
stagnated momentarily. In a specific embodiment, stagnation is achieved by closing a valve 
35 below the analyzer apparatus and diverting the flow stream around the sample chamber. 



30 



valid. 
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The time required for a measurement is almost entirely determined by the polarization 
time C'wait time*') of about 15 seconds. 

C. Hydrogen density measurements 

5 The hydrogen density or the total number of hydrogen atoms within tiie 

measurement volume is a by-product of any T| or T2 measurement It can be represented as 
the area tmder any Tj distribution and is typically normalized to the hydrogen density of a 
reference oil at measurement temperatures. At room temperature, the reference oil and 
water have the same hydrogen density. The reported mmiber is the relative hydrogen index 

10 Ih in the range 0-2, with accuracy aroimd 1%. 

Hydrogen density is automatically converted to hydrogen index Ih- The hydrogen 
index is hydrogen density relative to that of water at ambient conditions. Under the 
assumption that the oil contains only hydrogen and carbon atoms, the mass density p„^, the 
hydrogen index 1^, and the hydrogen-to-carbon ratio R are related as follows: 

15 

Ih - p„, 9R/(12+R) (6) 

See, for example, after Zhang et aL, 1998. Since the hydrogen index is measured, either the 
mass density or the H:C ratio can be computed jfrom an estimate of the other variable. 

20 It has been reported that most saturated hydrocarbon liquids have relative hydrogen 

indices of 1 within +/-5%. The hydrogen density in gases is significantly lower due to the 
overall lower density. Thus, a depressed hydrogen index serves as a first-order alert to the 
presence of gas and a change in the relationship between Tj and viscosity. Appel et aL 
(2000) reported a reduction of about 20% on live oil samples at 1 80°F. Under the 

25 assumption that all gas is methane (CH4), the observed hydrogen index can be approximated 
as follows: 

Ih = x(9/4p + (1-x)1, (6A) 

30 where x is the volumetric gas firaction (mVm^) and p, in g/cm^, is the density of methane. 
The density of methane follows from its temperature and pressure, and Eq. (6A) can be 
used to derive a first-order approximation for the gas firaction x. 



35 
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A Diffusion measurements 

Difiiision measurements in accordance with the present invention are implemented 
using steady-gradient spin-echo (SGSE) experiments. See, Kimmich et al^ 1997. The 
expejiment requires that the fluid flow be temporarily stopped. The idea of using the 

5 jftinges of a uniform-field volume for difiusometry is derived from so-called SSF-SGSE 
methods. Its main advantage over pulsed-field gradient spin-echo (PFGSE) difiRisometry is 
mstrumental simplicity and superior stability. The main drawback is a limit on sensitivity, 
which, for the downhole implementation, is approximately 10"^ cmVs. 

As noted above, the sensitive volume of the apparatus of this invention can be 

10 divided into an interior, homogeneous region and an exterior gradient region. The field in 
the fringe volume, which makes up about 1/3 of the total volume, can be approximated by a 
single field gradient value Gq. At short echo spacmgs (0.25 ms), the effect of the field 
gradient is too small to be relevant. The pulse sequence used both for diffusion 
measurements and for difiusivity calibration is shown in Fig. 7. 

15 In particular, two Carr-Purcell-Meiboom-Gill (CPMG) sequences with a short echo 

spacing (typically 0.25 ms) and a long spacing (TJ are alternated. The long echo spacing is 
selected as an integer multiple of the short spacing. In this case, echoes line up in time, i.e., 
occur at the same elapsed time since the excitation pulse and the ratio of their amplitudes 
can be formed. 

20 Assuming that the fluid under investigation has a Tj relaxation time, (the argument, 

also holds for an arbitrary distribution of Tj times) and diffiisivity D, the two echo trains for 
the short and the long echo spacing can be described as follows: 



25 



30 



35 



Ay = /^exp(-f/3^)and 

A, = I^K, exp(- t/T,) exp(- 1 / T^)+ 1^(1- K,) exp(- t/T,) 

where ' (7) 

l/T^^l/l2(yGJ,yD. 

The system parameter Ko is the gradient volume divided by the total volume. The 
hydrogen gyromagnetic ratio y is equal to 26,754 rad/s/gauss. Both Ko G© are 
temperature-dependent and are determined during calibration. The diffusivity D is derived 
from Eqs. 7 by taking the ratio of corresponding echoes, as follows: 

AJA, = K, exp(-f / T^) + (1- K^) (8) 
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This curve is fit to a uni-exponeatial model plus an offset In Fig, 8, the top two 
curves are the A, and signals for water at room temperature. Below is the ratio curve 
and the best-fit xmi-exponential model. Since D for water is known as 2.5x10'^ cmVs, these 
curves determine the calibration parameters Gq and Kq, * 

5 The top two curves in Fig. 8 are spin echo amplitudes at different echo spacings. 

The accelerated decay for the longer spacing is a manifestation of diffusion in the gradient 
region of the magnetic field. The ratio curve (below) is the sum of an exponential and a 
constant term, corresponding to the gradient-field region and the uniform-field region, 
respectively. The best-fit model curve is also plotted and is not distinguishable from the 

10 data. 

In a preferred embodiment, viscosity is determined as follows: 

;/ = 5xlO"^r/D (9) 

J J In this expression, the viscosity t] is measured in cp, the temperature T in Kelvin and 

the diffiisivity D in cmVs, The temperature may be obtained in a preferred embodiment 
from the RDT fluid temperature sensor. The proportionality factor has been determined by 
fitting Eq. (1) to data from pure alkanes and methane-alkane mixtures. 

2Q Calibration 

In a preferred embodiment, the system performs its own calibration, provided the 
sample chamber is filled with a known reference fluid (typically oil) and the system is 
heated through its operating temperature range. At specific temperature points, the bxiilt-in 
processor system records the sensor's resonance frequency and amplitude response and 

25 stores the results in a permanent calibration table. After calibration, the hydrogen dexisity is 
reported in percentage of that of the reference fluid, all times are in seconds and diffusivities 
are in cmVs. 



30 
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Applications 

Several practical applications are made possible by the apparatus and method of the 
present invention. Broadly, such applications relate to the measurement of relaxation 
distributions for the purpose of differentiation between reservoir hydrocarbons and oil- 

5 based mud (OBM) filtrates, as well as support for fluid-typing by wireline and LWD NMR 
tools. Other practical applications relate to calculation of hydrocarbon viscosity and gas-to- 
oil ratio and the modeling of flmd composition, as shown in the specific examples below. 

One suitable method for the identification of different oils is described in detail m 
U.S. Pat. No. 6,107,796, which is incorporated herein for all purposes. The differentiation 

10 is based on differences in distributions, where crude oils typically exhibit more complex 
distributions and characteristically short relaxation times. Similar differentiation is possible 
based on Tj relaxation measurements. It will be appreciated that in general relaxations 
are fiaster to measure and can be based on thousands of data pomts firom a CPMG sequence. 
Ti measurements are typically slower and yield few date points. However, as long as the 

15 underlying degrees of freedom in either distributions are not more than three-four, 

meaningful relaxation time distribution can be computed from about ten or more data points 
on a Tj recovery curve. Examples of specific practical applications follow. 

Connate oil v. mud filtrate differentiation 

20 Tj distributions can be used for qualitative fluid characterization of this type without 

invoking much of the relaxation theory. In a specific embodiment, the product tjTj is set to 
1 at T=300K and a simple viscosity scale is established from 1,000 cp to 1 cp, 
corresponding to the T, range 1 ms to 1,000 ms. See also Fig. 6. The water peak is 
distinctly offset due to the different molecular structure. The complexity of the distribution 

25 in the bottom panel of Fig. 6 is believed to reflect a distribution of internal mobilities due to 
a nux of short and long hydrocarbon chains. Contrast this appearance to the simple 
structure in the center panel in Fig. 6. This contrast provides a simple method to 
differentiate between complex oils (crudes) and simple oils (filtrates), that can be applied in 
a relatively straightforward manner using the techniques of the present invention. 

30 

NMR log interpretation 

Both wireline and LWD tools take their readings in the invaded zone, which is more 
or less flushed by filtrates. Porosity measurements by NMR are in fact hydrogen density 
readings in the fluid phase, calibrated to a water sample and corrected for temperature. 
35 Consequently, a lower-than-anticipated hydrogen index (dissolved gas, etc.) may cause an 
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undercall in porosity. Techniques are available to correct this effect for entrained methane, 
but it is more precise to sample the in^ed zone and directly determine the hydrogen index 
of the produced fluid. 

Fluid typing based on wireline/LWD NMR is complicated by the fact that T, and T2 

5 distributions are modulated by variations in pore sizes and by the bulk relaxation response 
of crade oils. The interpretation makes simplifying assumptions about the hydrocarbon 
phase (non-wetting, single T^ for oil; single T, for gas). In accordance with the present 
invention this approach can be refined by determining the actual relaxation profiles fi:om 
fluid samples and by feeding this information back into the saturation calculation. 

1 0 Wireline and LWD NMR tools operate at frequencies between 0.5 MHz and 2 MHz, 

while the fluid analyzer described above operates at about 4 MHz. It will be appreciated 
that this difference is irrelevant for the relaxometry of fluids with Ti s Tj. This condition 
indicates a uniform and flat energy spectrum in the employed firequency range. 

15 Hydrocarbon viscosity 

Determination of viscosity under true reservoir conditions has always been a 
challenge. As noted above, the NMR-derived diffiisivity D has a universal correlation with 
viscosity. In order to obtain a continuous viscosity reading, however, it is desirable to also 
derive a viscosity estimate from Tj and 1^ alone. For the case of no dissolved hydrocarbon 

20 gas (i.e. close to 1) and no entrained paramagnetic oxygen, currently the best available 
correlation is: 



25 where T, and Tj are in s, T is in K, and t] is in cp (See, Lo et a/., 2000). To compute 
viscosity from a T| distribution, one can use the approximation: 



30 where T,^"* is the geometric mean taken over the T, distribution. It is expected that this 
correlation can be improved by includmg the hydrogen index once a sufficient body of 
downhole NMR data becomes available. 



Tt -T2 « 9.6 10-^ 



T/tj (degassed alkanes), 



(11) 



Ti« 9.6 10-^T/Ti^"' 



(12) 



Gas/oil ratio (GOR) 



35 



-23- 



wo 02/079808 



PCTAJS02/09819 



For dead oils (GOR=0), the relationship between and D is basically linear, as can 
be seen with reference to Eqs. (1) and (2). Increasing gas contents introduces a deviation 
fit>m the linear behavior. Modification of Eq. 2 and substituting from Eq. 1 leads to: 



where f(GOR) is a function tiiat has been detemuned empirically for methane-n-alkane 
mixtures. See Freedman et al., 2000: The proportionality factor in Eq. (13) is 
approximately 2x10^ (Lo et al, 2000). Accordingly, the GOR can be determined 
10 graphically by means of a crossplot Tj v. D with GOR as parameter, as shown in Fig. 1 3 of 
the Lo et al reference.. 

Emerging practical applications aim at a comprehensive imderstanding of the 
reservoir fluid system, the evaluation of which can be improved using the present invention. 

1 5 Although the present invention has been described in connection with the preferred 

embodiments, it is not intended to be limited to the specific form set forth herein, but on tiie 
contrary, it is intended to cover such modifications, alternatives, and equivalents as can be 
reasonably included within the spirit and scope of the invention as defined by the following 
claims. 



5 



T„T2ccD/f(G0R), 



(13) 
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WHAT IS CLAIMED IS! 

1 . A method for analyzing fonnation fluids in a borehole environment, 
comprising the steps of: 

(a) introducing formation fluids in a flow-through vessel located m the borehole; 
5 (b) generating a substantially uniform static magnetic field in the vessel with a 

defined magnetic field direction; 

(c) generating pulsed oscillating magnetic fields for exciting nuclei of formation 
fluids in a first portion of the vessel, said oscillating fields having a magnetic direction 
substantially perpendicular to the direction of the static magnetic field; 
1 0 (d) receiving nuclear magnetic resonance (NMR) relaxation signals from excited 

nuclei of the fluid in a second portion of the vessel, smaller than the first portion; and 

(e) analyzing the received signals to determine fluid properties at any flow rate 
of the formation fluids in the vessel below a predetermined non-zero threshold. 

1 5 2. The method of claim 1 , wherein in step (c) the pulsed magnetic fields are 

according to a saturation recovery pulse sequence. 

3. The method of claim 2, wherein an excitation pulse in step (c) excites nuclei 
in substantially the entire first portion of the vessel. 

20 

4. The method of claim 1, wherein the threshold value for the flow rate is 
determined by the length of the first portion. 

5. The method of claim 1 , wherein the second portion of the vessel is 
25 positioned downstream from the first portion of the vessel. 

6. The method of claim 5, wherein the first and second portions of the vessel 
overlap at least in part. 

30 7. The method of claim 1 , wherein steps (a)-(d) are performed substantially 

continuously. 

8. The method of claim 7, wherein the step of analyzing comprises assessing 
mud filtrate contamination for the formation fluids. 
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9. The method of claun 7, v^erem the step of analyzing comprises monitoring 
the T, pro£U[e of fluids passed through the vessel. 

1 0. The method of claim 1 , wherein parameters of the pulsed oscillating 

S magnetic fields are selected to effect saturation in the first and second portions of the vessel. 

11. A method for analyzing fluids, comprising the steps of: 

(a) providing a flow-through passage for fluids in a measurement chamber, 
the fluids having flow rate within a pre-detemiined range; 
10 (b) performing a nuclear magnetic resonance (NMR) experiment to excite 

substantially all nuclei of the fluids in the chamber at a given time interval; and 

(c) processing NMR signals obtained from a portion of the chamber in said 
experiment without regard for the flow rate of the fluids in the chamber. 

1 5 12. The method of claim 1 1 , wherein the upper limit of the range in step (a) is 

determined based on the length of the chamber. 

13 . The method of claim 1 1 , wherein the flow-through passage in step (a) has an 
inlet end and an outlet end and at least a portion of the passage between the inlet and outlet 

20 ends is widened. 

14. The method of claim 1 1 further comprising, prior to step (b), the step of 
subjecting fluids in the flow-through passage to a prepolarization static magnetic field to 
polarize hydrogen nuclei. 

25 

15. The method of claim 1 1, wherein the prepolarization magnetic field is 
stronger tiian static magnetic fields used in the NMR experiment in step (b). 

1 6. The method of cledm 1 1 , wherein the NMR experiment in step (b) is 
30 performed according to a saturation recovery pulse sequence. 

1 7. The method of claim 1 1, wherein an excitation pulse of the saturation 
recovery pulse sequence excites nuclei in substantially the entire measurement chamber. 
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1 8. The method of claim 1 1 , wherein the step of processing comprises 
monitoring the T{ profile of fliiids passed through the measurement chamber. 

19. The method of claim 11, wherein steps (a)-(c) are performed substantially 
S continuously. 

20. The method of claim 1 9, wherem the fluids are geological formation fluids m 
a borehole, and the step of processing comprises assessing mud filtrate contamination for 
the formation fluids. 

10 

21 . An apparatus for analyzing formation fluids in a borehole environment, 
comprising; 

(a) means for introducing formation fluids in a flow-through vessel located in 
the borehole; 

15 (b) means for generating a substantially uniform static magnetic field in the 

vessel with a defined magnetization direction; 

(c) means for generating pulsed oscillating magnetic fields for exciting nuclei of 

formation fluids in a first portion of the vessel, said oscillating fields having a magnetic 

direction substantially perpendicular to the durection of the static magnetic field; 
20 (d) means for receiving nuclear magnetic resonance (NMR) relaxation signals 

from excited nuclei of the fiuid in a second portion of the vessel, smaller than the first 

portion; and 

(e) means for analyzing the received signals to determine fluid properties at any 
flow rate of the formation fluids in the vessel below a predetermined non-zero threshold. 

25 

22. The apparatus of claim 21, wherein the means for generating pulsed 
oscillating magnetic fields comprises means for generating a saturation recovery pulse 
sequence. 

30 23. The apparatus of claim 21, wherein the second portion of the vessel is 

positioned downstream from the first portion of the vessel. 

24. Hie apparatus of claim 23, wherein the first and second portions of the vessel 
overlap at least in part. 
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25. An apparatus for analyzing fluids, comprising: 

(a) means for providing a flow-through passage for fluids in a measurement 
chamber, the fluids having flow rate within a pre-detennined range; 

(b) means for performing a NMR experiment to excite substantially all 
5 nuclei of the fluids in the chamber at a given time interval; and 

(c) means for processing NMR signals obtained fix)m a portion of the 
chamber in said experiment without regard for the flow rate of the fluids in the chamber. 

26. The apparatus of claim 25, wherein the flow-through passage in (a) has an 

1 0 inlet end and an outlet end and at least a portion of the passage between the inlet and outlet 
ends is widened. 

27. The apparatus of claim 25 further comprising means for subjecting fluids in 
the flow-through passage to a prepolarization static magnetic field. 

15 

28. The apparatus of claim 27, wherein the prepolarization means comprises a 
permanent magnet assembly enclosing the passage and having a magnetic field stronger 
than static magnetic fields generated by the means for perforaiing NMR experiment. 

29. An apparatus for analyzing downhole formation fluids, comprising: 

(a) a conduit for introducing formation fluids into the apparatus and for 
providing flow-through passage, the conduit having an inlet end and an outlet end; 

(b) at least one magnet assembly enclosing the conduit for generating in the 
conduit a substantially imiform static magnetic field with a defined magnetic field direction; 

(c) at least one transmitting antenna operative to generate pulsed magnetic fields 
in the conduit in a direction substantially perpendicular to the static field direction for 
exciting nuclei of fluids contained in the conduit; and 

(d) at least one receiving antenna operative to receive NMR signals firom fluids 
in a portion of the conduit, the receiving anterma being shorter than the at least one 
transnaitting coil, so that received NMR signals correspond only to a portion of the excited 
nuclei. 

30. The apparatus of claim 29, wherein the portion of the conduit between the 
inlet and outlet ends is wider in dimension than the conduit at either end. 
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3 1 . The apparatus of claim 3 0, M^ierein the conduit has cylindrical shape and the 
diameter of the wider portion of the conduit is approximately 2 cm, and the diameter of the 
inlet and outlet ends is approximately O.S cm. 

5 32. The apparatus of claim 30 further comprising a diffuser positioned near the 

inlet end of the conduit, for providing consistent fluid flow velocity over the wider portion 
of the conduit. 

33. The apparatus of claim 29, wherein the conduit comprises thin collimator 
1 0 tubes, for slowing down fluid entering the inlet end of the conduit. 

34. The apparatus of claim 29, wherein the conduit is adapted for attachment to 
the flow line of a modular wireline logging tool. 

15 35 . The apparatus of claim 29, wherein the conduit is made of one or more of: 

ceramic, glass or PEEK material. 

36. The apparatus of claim 35, wherein the conduit is enclosed in a Faraday 
shield slotted along the length of the conduit. 

37. The apparatus of claim 29, wherein the at least one magnet assembly has a 
polarization portion located near the inlet end, and a resonance portion located near the 
outlet end, the strength of the magnetic field in the polarization portion being higher than 
the strength of the magnetic field in the resonance portion. 

38. The apparatus of claim 37,' wherein the at least one magnet assembly further 
comprises a stabilization portion positioned between the polarization and the resonance 
portions, for allowing magnetic spins in the fluid flow following the polarization portion to 
settle to an equilibrium polarization. 

39. The apparatus of claim 37, wherein the polarization portion of the magnet 
assembly has a Halbach configuration. 
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40. The apparatus of claim 39, the Halbach magnet of the polarization 

portion is dimensioned to produce a magnetic field of approxunately 2,000-2,500 Gauss 
within the conduit. 

5 41 . The apparatus of claim 3.9, wherein the Halbach magnet of the polarization 

portion is made of low-tempemture coefQcient Samarium-Cobalt material. 

42. The apparatus of claim 38, wherein the stabilization portion of the magnetic 
assembly comprises a plurality of magnetic segments, which are magnetized in the same 
10 direction. 



43 . The apparatus of claim 29, wherein the at least one transmitting anteima is a 
solenoidal coil wound around a measurement portion of the conduit. 

1 5 44. The apparatus of claim 43 further comprising a shield moxmted between the 

transmitting antenna and a portion of the magnet assembly. 

45. The apparatus of claim 29, wherein application of a jfrequency-swept pulse in 
the transmitting antenna saturates nuclear magnetization in a predetermined portion of the 

20 conduit. 

46. The apparatus of claun 29, wherein at least one receiving anteima is a 
solenoidal coil wound around a portion of the magnet assembly. 

25 47. The apparatus of claim 46, wherein a portion of the magnet assembly 

corresponding to the receiving antenna has magnetic field variations less than about 0.1%. 

48. The apparatus of claim 29, wherein at least one transmitting anteima is a 
solenoidal coil, at least one receiving anteima is a solenoidal coil, and the receiving coil is 

30 positioned closer to the outiet end of the conduit 

49. The apparatus of claim 29, wherein the at least one transmitting and at least 
one receiving antennas overlap at least in part, and the receiving antenna is positioned closer 
to the outiet end of tiie conduit. 
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50. The apparatus of claim 49, wherein the transmitting antenna is 2 to 3 times 
longer than tiie receiving antenna. 

5 1 . The apparatus of claim 48, further comprising a temperature compensating 
5 circuit attached to the at least one transmitting and at least one receiving antenna coils. 

52. A tester module for use with modular downhole formation testers for 
downhole nuclear magnetic resonance (NMR) testing of formation fluids comprising: 

a vessel for providing flow-through passage of formation fluids and for 
10 conducting downhole NMR measurements, said vessel being adapted to withstand borehole 
environment conditions; 

at least one tubular magnet defining a longitudinal axis, the magnet having 
magnet sections with magnetization direction(s). perpendicular to the longitudinal axis, and 
enclosing the passage to generate therein a static magnetic field with predeteraiined 
15 magnetic field direction; 

at least one radio firequency (RF) transmitter operative to generate pulsed RF 
magnetic fields in a first portion of the passage in a direction substantially perpendicular to 
the static field direction for exciting nuclei of fluids in the passage; 

a receiver for acquiring NMR signals from excited nuclei in a second portion 
20 of the passage smaller than the first portion and positioned downstream therefirom, and 

a processor for analyzing properties of fluids in the passage based on signals 
from the receiver. 
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